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proposed for light propagation in fibre Bragg gratings (FBGs) with phase defects. A periodic sequence of 
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mapped into the spectral response of the FBG. For the semi-infinite FBG Tamm surface states can appear 
and can be visualized as narrow resonance peaks in the transmission spectrum of the grating. 
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1. Introduction 

The Kronig-Penney model for the non-relativistic Schrodinger equation [1] is a well-known model in solid-state 
physics that describes the electronic band structure of an idealized one-dimensional crystal. The Kronig-Penney 
model has served on many occasions as a paradigmatic model to study a wide variety of physical phenomena, 
including band structure properties, localization effects in disordered lattices, electronic properties of superlattices, 
Peierls transitions and quark tunnelling in one-dimensional nuclear models. Relativistic extensions of the Kronig- 
Penney model (also referred to as the Dirac-Kronig-Penney model) have been discussed by several authors (see, 
for instance, [2-13] and references therein), and the impact of relativity on the band structure and localization, 
such as shrinkage of the bulk bands with increasing band number, have been highlighted on many occasions. In 
earlier studies, the Dirac-Kronig-Penney model also attracted some attention and caused a lively debate about 
the existence of so-called Dirac surface states, i.e. relativistic surface Tamm states which disappear in the non- 
relativistic limit [3, 7, 14-19]. In recent years, there has been an increased interest in simulations of relativistic 
quantum effects using different physical set-ups, and analogues of such fundamental phenomena as Zitterbewegung 



* E-mail: longhi@fisi.polimi.it 



Photonic realization of the rclativistic Kronig- Penney model and relativistic Tamm surface states 



and Klein tunnelling -rooted in the Dirac equation- have been proposed and demonstrated for electrons in graphene 
[20, 21] and for matter waves using trapped cold atoms [22]. Light propagation in guiding optical structures has 
been also shown to provide a beautiful laboratory system to investigate the classical analogous of a wide variety 
of coherent non-relativistic [23] and relativistic [24-28] quantum phenomena. In particular, it was recently shown 
that light propagation in fibre Bragg gratings (FBGs), i.e. optical fibres with a superimposed modulation of the 
refractive index profile, provides an experimentally accessible laboratory tool to simulate in optics the massive 
one-dimensional Dirac equation, and a photonic realization of the Dirac oscillator (i.e. the relativistic extension 
of the quantum harmonic oscillator) has been proposed [29]. 

In this work we propose a photonic realization of the relativistic Kronig-Penney model and relativistic surface 
Tamm states based on light propagation in FBGs with phase defects. Light propagation in a FBG with a periodic 
sequence of phase slips is shown to simulate the relativistic Kronig-Penney model, the relativistic band structure 
of which being mapped into the spectral transmission of the FBG. Similarly, a semi-infinite FBG with phase 
defects interfaced with a uniform FBG with a different modulation period is shown to support Tamm surface 
states analogous to the relativistic Tamm states. Such surface states are responsible for narrow resonance peaks in 
the transmission spectrum of the grating. The paper is organized as follows. In Section 2, the photonic realization 
of the Dirac-Kronig-Penney model for an infinitely-extended crystal, based on a superstructure FBG, is presented. 
The photonic analogues of relativistic Tamm states in a semi-infinite lattice model are discussed in Section 3, and 
shown to appear as narrow resonance peaks in the spectral transmission of the FBG. Finally, the main conclusions 
are outlined in Section 4. 

2. Photonic realization of the Dirac-Kronig-Penney model 

In this section a photonic realization of the Dirac-Kronig-Penney model for an infinitely-extended lattice, based 
on Bragg scattering of light waves in a superstructure FBG comprising a periodic sequence of phase slips, is 
proposed. Let us consider light propagation in a FBG with a longitudinal effective refractive index given by 



where no is the effective mode index in absence of the grating, An <^ no is the peak index change of the grating, 
A is the nominal period of the grating defining the reference frequency ujb = 7rc/(Ano) of Bragg scattering, c is 
the speed of light in vacuum, and m{z), (f>[z) are the slow variation, as compared to the scale of A, of normalized 
amplitude and phase, respectively, of the index modulation. The periodic index modulation of the grating leads 
to Bragg scattering between two counterpropagating waves at frequencies close to ojs. By letting 



n{z) — no + An m{z) cos 



2nz 



(1) 



E{z, t) = t) exp [—iujBt + iksz + i(j}{z)/2] + ip2{z, t) exp [—iujBi — iksz — i(j}{z)/2] + c.c. 



(2) 




Stefano Longhi 



for the longitudinal electric field amplitude in the fibre, where fcs — n/A, the slowly-varying envelopes ipi and if2 
of counterpropagating waves satisfy the coupled-mode equations (see, for instance, [30]) 
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(3) 
(4) 

(5) 



and Vg ~ c/na is the group velocity at the Bragg wavelength \b = 2-nc/u)B- The analogy between Bragg scattering 
of counterpropagating light waves in the FBG and the temporal evolution of two-component spinor Dirac equation 
in presence of an electrostatic potential is at best captured by introducing the dimensionless variables x = z/Z 
and T = t/T , with characteristic spatial and time scales defined by 



2no 
fcs An 



r = 



After introduction of new envelopes 

^1,2(2:) = 

Eqs.(3) and (4) can be cast in the Dirac form 



V2 



.dip . dip , , , , , , 

= —iai—ip + m{x)a3ip + V(x)ip 
OT ax 



(6) 



(7) 



(8) 



for the spinor wave function ip = {ipi,ip2)'^, where 



V{x)^ 



2 dx 



(9) 



and (Ti,3 are the Pauli matrices, defined by 



0"! = 



1 

1 



(T3 = 



1 

-1 



(10) 



In its present form, Eq.(8) is analogous to the one-dimensional Dirac equation, written in natural units ?i = c = 1, 
in the presence of an external electrostatic potential V{x), m{x) playing the role of a dimensionless (and generally 
space-dependent) rest mass. 

The Dirac-Kronig-Penney model for an infinitely-extended lattice corresponds to a constant mass ■m{x) = mo and 
to a potential V{x) given by the superposition of equally-spaced J-like barriers, namely [2] 



V{x) = Vo J2 " 



(11) 
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where Vb > is the area of the barrier and a is the lattice period. Stationary solutions i>{x, t) — tl^oix) exp(— iiSr) 
to the Dirac equation (8) in the periodic potential (11) with energy E are of Bloch-Floquet type, i.e. tf}o{x + a) — 
ipo{x) exp{iqa), where q is the Bloch wave number which varies in the first Brilloiun zone [—it /a < q < n/a). The 
corresponding energy spectrum is composed by a set of allowed energy bands E — E{q), which are defined by the 
following implicit equation (see, for instance, [17]) 

cos(ga) = cos(Vb) cos(fi:a) H sin(Vb) sin(Ka), (12) 

where we have set 

K = ^E'^ -ml. (13) 

Equation (12) defines the dispersion relation of the relativistic Kronig-Penney model, which has been investigated 
by several authors (see, for instance, [2, 4, 7]). The ordinary non-relativistic limit of the Kronig-Penney model 
is attained from Eqs.(12) and (13) for Vb ^ 1 and for energies E close the mo, for which the energy-momentum 
relation (13) reduces to the non-relativistic one [E ~ mo -I- K^/(2mo)]; in this regime, the dispersion relation (12) 
reduces to cos(ga) = cos(Ka) 4- (moVb/ft;) sin(Ka), which is the ordinary dispersion relation encountered in the 
non-relativistic Kronig-Penney model. For larger energies E but still for a low barrier area Vo <^1, non-relativistic 
effects come into play as perturbative effects, which modify positions and widths of the allowed energy bands. 
Non-relativistic effects deeply modify the band structure of the crystal for potential strengths Vb of the order ~ 1. 
In particular, if Vb is an integer multiple of tt, all band gaps disappear and the dispersion relation reduces to the 
one of a relativistic free particle [Eq.(8) with V{x) = Q], as if the 5-barriers were absent. 

In our photonic context, owing to Eq.(9) the Dirac-Kronig-Penney model simply corresponds to an infinitely- 
extended uniform FBG with a superimposed periodic sequence of lumped phase slips of equal amplitude = 2Vb 
and spaced by the distance a. The circumstance that the effects of the 5 barriers disappear in the Dirac-Kronig- 
Penney model when Vb is an integer multiple tt is simply due to the fact that, under such a condition, the phase 
slips are integer multiplies of 2-n, and thus the grating has no phase defects and mimics the dynamics of a one- 
dimensional free relativistic Dirac particle. It should be noticed that, in another physical context, superstuctured 
FBGs comprising a periodic sequence of tt phase slips, corresponding to the special case Vb = 7r/2, have been 
recently proposed and demonstrated to realize light slowing down [31, 32]; however, their connection to the Dirac- 
Kronig-Penney model was not noticed. It should be also noticed that any FBG structure has a finite length 
L, and thus some kind of lattice truncation should be introduced in the idealized model. As shown e.g. in 
[31, 32], apodization of the FBG amplitude profile m(a;), obtained by slowly decreasing m{x) from its constant 
value m{x) — mo in the uniform grating region to zero at x — >■ ±oo, enables to adiabatically inject and eject light 
waves in the uniform grating region avoiding truncation effects. Correspondingly, the band structure of the Dirac- 
Kronig-Penney model is mapped into the alternation of stop/transmission bands observed in spectrally-resolved 
transmission measurements of the FBG. The spectral transmission (t) and reflection (r) coefficients for the Dirac 
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equation with m{x) — > at a; — >■ ±00, corresponding to left wave incidence, are defined from the scattering 
solution to (8) with the following asymptotic behavior 



r / 
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4 



exp[— i_Ea; — iEt + i(f>{x)/2] 



X — >■ —00 



V 



exp[iEx — iEr — i(fi{x) /2] 



(14) 



Note that, in our photonic analogue, the energy E of the scattered waves corresponds to the (normalized) frequency 
detuning of light waves, propagating in the grating, from the Bragg reference frequency lub- 

As an example. Fig. 1(c) shows the numerically-computed spectral power transmission \t{E)\^ (in dB units) for a 
uniform FBG comprising a sequence of tt phase slips [see Fig. 1(b)] and with a super-Gaussian apodizion profile 
m{x) shown in Fig. 1(a); parameter values are Vo ~ 7r/2, a = 2, mo = 1 and L — 50. The power spectral 
transmission has been computed using a standard transfer matrix method (see, for instance, [30]). The shaded 
areas in Fig. 1(c) indicate the stop bands of the ideal (infinitely-long) Dirac-Kronig-Penney lattice as predicted 
by Eq.(12). In physical units, for typical parameter values no = 1.45, An = 1 x 10~*, and As — 1560 nm, which 
are consistent with the those of Ref.[32] , the spatial Z and temporal T scales for the example in Fig.l are Z ~ 5 
mm and T ~ 24 ps, respectively. Hence, in physical units the grating length is L ~ 25 cm, the phase slips are 
separated by the distance aZ ~ 1 cm, and the unit scale of (nonangular) frequency detuning E from the reference 
Bragg frequency in Fig. 1(c) is l//{2nT) ~ 6.6 GHz. 

3. Relativist ic Tamm surface states 



The existence of surface Tamm states for the relativistic Kronig-Penney model has been widely studied in earlier 
papers by several authors [3, 7, 14-19]. In such studies, a debate was raised about the proper boundary conditions 
that should be imposed to the relativistic wave function at a (5 barrier (see, for instance, [17, 34, 35]). As earlier 
works [3, 14-16], based on incorrect boundary conditions, suggested that the relativistic treatment yields a new 
class of surface states (the so-called Dirac surface states) which do not correspond the common Tamm states [33] 
in the non-relativistic limit, it was subsequently realized that application of more physical boundary conditions 
does not yield any surface state which violates the Tamm condition in the non-relativistic limit [17]. In our 
photonic system, the appropriate boundary conditions in presence of discontinuities of the phase (j>{^) [sjad hence 
a (5-like behavior of the potential V{x); see Eq.(9)] are readily obtained from Eqs.(3) and (4), and coincide with 
those used by Subramanian and Bhagwat in the study of relativistic Tamm states [17]. The relativistic extension 
of the Tamm model is defined by the potential (see, for instance, [7, 17]) 



V{x) = . 



Vi X <0 

VoJ2^=i S(x - ria) x>0 



(15) 
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Figure 1. Photonic analogue of the Dirac-Kronig-Pcnney model in a superstructure FBG comprising a periodic sequence of phase 
slips, (a) and (b): amplitude and phase profiles of the grating, (c) Numerically-computed spectral power transmission. 
Parameter values are given in the text. The dashed areas in (c) are the stop bands of the corresponding Dirac-Kronig- 
Pcnney infinite lattice. 



Surface states are found as localized solutions to Eq.(8), near the surface a: = 0, satisfying the appropriate 
boundary conditions. The energies of such states, if any, should obviously fall in a gap of the energy spectrum of 
the infinitely-extended Dirac-Kronig-Penney model. As in Ref. [17], we limit to consider the case Vi < mo, for 
which surface states occur at the energies E in the interval (mo, mo + Vi) satisfying the equation 



/^cotg(/^a) = ViCotg(Vb) — K 



(16) 



provided that E — Vi — i^cotg(V()) > [17], where k is defined by Eq.(13) and 



K 



(17) 



In our photonic system, the potential V{x) defined by Eq.(15) and supporting the surface Tamm states at the 
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X = boundary corresponds to a phase profile ^(a;) of tfie FBG wliicli is composed by a linearly increasing 
part with a slope 2Vi for a; < 0, and by a staircase phase profile for a; > [see, for instance, Fig. 2(b) to be 
commented below], as one can see after integration of Eq.(9), <j){^) = 2 d^V{^). Physically, such a FBG 
basically corresponds to two adjacent sections of uniform grating regions but with different grating periods, with 
the second grating region (at a; > 0) comprising a sequence of equally-spaced phase slips, at a distance a, equal 
to Aff) — 2Vo- Similarly to the photonic realization of the Dirac-Kronig-Penney model discussed in the previous 
section, a finite grating length is introduced by apodization of the amplitude profile m(x). In this case, the surface 
states actually become resonances of the FBG, i.e. trapped light states near x — which are weakly coupled 
to the external regions of the fibre because of evanescent tunnelling (see, for instance, [29]). The existence of 
surface states can be thus simply recognized by the appearance of narrow resonance peaks embedded in a stop 
band region of the transmission spectrum of the grating. An example of the photonic analogue of relativistic 
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Tamm states in a FBG is shown in Fig. 2, corresponding to parameter values Vb = 7r/2, a = 2, mo = 1, (same 
as in Fig.l) and Vi = 0.8. For such parameter values, a numerical analysis of the Tamm condition [Eq.(16)] 
indicates that there exists one Tamm surface state at the energy value E = Eo — 1.474, and thus a resonance 
peak in the spectral power transmission of the grating is expected at such a detuning value. Figures 2(a) and 
2(b) show the amplitude [Fig. 2(a)] and phase [Fig. 2(b)] profiles of the FBG, whereas the corresponding power 
transmission (numerically-computed by the transfer matrix method) is depicted in Fig. 2(c). Note that, according 
to the theoretical analysis, a strong and narrow resonance peak, at the detuning E — Eq, is clearly observable in 
the transmission spectrum [see the inset of Fig. 2(c)], which is the signature of the surface state localized between 
the two grating regions at around x — 0. For parameter values no = 1.45, An = 1 x 10^*, and As — 1560 nm (as 
in Fig.l), the spatial Z and temporal T scales are given by Z ~ 5 mm and T ~ 24 ps, respectively, and hence in 
physical units the grating length is L ~ 25 cm, the phase slips are separated by the distance aZ ~ 1 cm, and the 
unit scale of (nonangular) frequency detuning E from the reference Bragg frequency in Fig. 2(c) is 1// {2nT) ~ 6.6 
GHz. 

4. Conclusions 

In conclusion, in this work a photonic realization of the Dirac-Kronig-Penney model, describing the band 
structure of a periodic potential in the relativistic regime, as well as of relativistic surface Tamm states, 
has been proposed. Our photonic analogue is based on Bragg scattering of light waves in a uniform FBG 
with a periodic sequence of phase slips. Band structure and surface states of the relativistic lattice model 
can be simply observed from spectrally-resolved transmission measurements of the FBG. Design parameters 
of the grating structures, which are compatible with the current FBG writing technology [32], have been presented. 

This work was supported by the Italian MIUR (PRIN-2008 project "Analogic ottico-quantistiche in strutture 
fotoniche a guida d'onda"). 

References 



[1] R. de L. Kronig and W.G. Penney, Proc. Roy. Soc. London Ser. A 130, 499 (1931). 

[2] F. Dominguez-Adame, Am. J. Phys. 55, 1003 (1987). 

[3] S.G. Davison and M. Streshcka, J. Phys. C 2, 1802 (1969). 

[4] B. H. J. McKellar and G.J. Stephenson, Phys. Rev. A 36, 2566 (1987). 

[5] F. Dominguez-Adame, J. Phys.: Condes. Matter 1, 109 (1989). 

[6] I.M. Mladenov, Phys. Lett. A 131, 313 (1989). 

[7] G.J. Clerck and B.H.J. McKellar, Phys. Rev. C 41, 1198 (1990). 



Stefano Longhi 



[8] C.L. Roy and C. Basu, J. Phys. Chem. Solids 52, 745 (1991). 

[9] F. Dominguez-Adame and A. Sanchez, Phys. Lett. A 159, 153 (1991). 

[10] G.J. Clerck and B.H.J. McKellar, Phys. Rev. B 47, 6942 (1993). 

[11] C. Basu, C.L. Roy, E. Macia, F. Dominguez-Adame, and A. Sanchez, J. Phys. A 27, 3285 (1994). 

[12] F. Dominguez-Adame, E. Macia, A. Khan, and C.L. Roy, Physica B 212, 67 (1995). 

[13] M. Barbier, F.M. Peeters, P. Vasilopoulos, and J.M. Pereira, Phys. Rev. B 77, 115446 (2008). 

[14] M.L. Glasser and S.G. Davison, Int. J. Quantum Chem. 4, 867 (1970). 

[15] S.G. Davison and J.D. Levine, Solid State Phys. 25, 32 (1970). 

[16] M. Steslicka and S.G. Davison, Phys. Rev. B 1, 1858 (1970). 

[17] R Subramanian and K.V. Bhagwat, J. Phys. C 5, 798 (1972). 

[18] S. Yuyi, Surface Science 108, L477 (1981). 

[19] C.L. Roy and J.S. Pandey, Physica 137A, 389 (1986). 

[20] A.F. Young and P. Kim, Nature Phys. 5, 222 (2009). 

[21] N. Slander, B. Huard, and D. Goldhaber-Gordon, Phys. Rev. Lett. 102, 026807 (2009). 

[22] R. Gerritsma, G. Kirchmair, F. Zhringer, E. Solano, R. Blatt, and C. F. Roos, Nature 463, 68 (2010). 

[23] S. Longhi, Laser Photonics Rev. 3, 243 (2009). 

[24] R. A. Sepkhanov, Ya. B. Bazaliy, and C. W. J. Beenakker, Phys. Rev. A 75, 063813 (2007). 

[25] X. Zhang, Phys. Rev. Lett. 100, 113903 (2008). 

[26] L.-G. Wang, Z.-G. Wang, and S.-Y. Zhu, EPL 86, 47008 (2009). 

[27] S. Longhi, Phys. Rev. B 81, 075102 (2010). 

[28] O. Bahat-Treidel, O. Peleg, M. Grobman, N. Shapira, T. Pereg-Barnea, and M. Segev, Phys. Rev. Lett. 104, 

063901 (2010). 

[29] S. Longhi, Opt. Lett. 35, 1302 (2010). 

[30] T. Erdogan, J. Lightwave Technol. 15, 1277 (1997). 

[31] D. Janner, G. Galzerano, G. Delia Valle, P. Laporta, S. Longhi, and M. Belmonte, Phys. Rev. E 72, 056605 
(2005). 

[32] S. Longhi, D. Janner, G. Galzerano, G. Delia Valle, D. Gatti, and P. Laporta, Electron. Lett. 41, 1075 (2005). 

[33] I. Tamm, Z. Phys. 76, 849 (1932). 

[34] B. Sutherland and D.C. Mattis, Phys. Rev. A 24, 1194 (1981). 

[35] C.L. Roy, Phys. Rev. A 74, 3417 (1993). 



9 



